Splitting of a Cooper pair by a pair of Majorana bound states 
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Majorana bound states are spatially localized superpositions of electron and hole excitations in 
the middle of a superconducting energy gap. A single qubit can be encoded nonlocally in a pair 
of spatially separated Majorana bound states. Such Majorana qubits are in demand as building 
blocks of a topological quantum computer, but direct experimental tests of the nonlocality remain 
elusive. Here we propose a method to probe the nonlocality by means of crossed Andreev reflection, 
which is the injection of an electron into one bound state followed by the emission of a hole by the 
other bound state (equivalent to the splitting of a Cooper pair over the two states). We have found 
that, at sufficiently low excitation energies, this nonlocal scattering process dominates over local 
Andreev reflection involving a single bound state. As a consequence, the low-temperature and low- 
frequency fluctuations Sli of currents into the two bound states i = 1, 2 are maximally correlated: 
ShSh = Iff. 

PACS numbers: 03.75.Lm, 73.21.-b, 74.45. +c, 74.78.Na 
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Majorana bound states are coherent superpositions of 
electron and hole excitations of zero energy, trapped 
in the middle of the superconducting energy gap by 
a nonuniformity in the pair potential. Two Majorana 
bound states nonlocally encode a single qubit (see Fig. 
[3 top panel). If the bound states are widely separated, 
the qubit is robust against local sources of decoherence 
and provides a building block for topological quantum 
computation [3, 0. 

While Majorana bound states have not yet been 
demonstrated experimentally, there is now a variety of 
candidate systems. In an s-wave superconductor, zero- 
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FIG. 1: Top panel: Energy diagram of two Majorana bound 
states (levels at zero energy), which split into a pair of levels 
at ±Em upon coupling. Whether the upper level is excited 
or not determines the states |1) and |0) of a qubit. Crossed 
Andreev reflection probes the nonlocality of this Majorana 
qubit. Lower panel: Detection of crossed Andreev reflection 
by correlating the currents Li and 1% that flow into a super- 
conductor via two Majorana bound states. 



point motion prevents the formation of bound states at 
zero energy. Early proposals for Majorana bound states 
therefore considered p-wave superconductors [3|, \A\ , with 
Sr2Ru04 as a candidate material @, or p-wave superflu- 
ids formed by fermionic cold atoms Q. More recently, it 
was discovered @, H, [1] that Majorana bound states can 
be induced by s-wave superconductivity in a metal with 
a Dirac spectrum (such as graphene or the boundary of 
a topological insulator). Several tunneling experiments 
have been proposed [l(J, EH E3| to search for the Majo- 
rana bound states predicted to occur in these systems. 

Here we show that crossed Andreev reflection [rJ EH 
EE[ by a pair of Majorana bound states is a direct probe 
of the nonlocality. Crossed Andreev reflection is the non- 
local conversion of an electron excitation into a hole ex- 
citation, each in a separate lead. Local Andreev reflec- 
tion, in contrast, converts an electron into a hole in the 
same lead. Equivalently, local Andreev reflection injects 
a Cooper pair in a single lead, while crossed Andreev 
reflection splits a Cooper pair over two leads. We have 
found that at sufficiently low excitation energies, local 
Andreev reflection by a pair of Majorana bound states is 
fully suppressed in favor of crossed Andreev reflection. 

The suppression is not a property of the dispersion 
relation in the leads (as in Refs. [lg, E3)j but directly 
probes the Majorana character of the Hamiltonian Q 



Hm = i-EA/7172, 



(1) 



of the pair of weakly coupled bound states (labeled 1 and 
2). The Ti's are Majorana operators, defined by 7; = 7J, 
lilj+ljli = Zfiij- The coupling energy Em splits the two 
zero-energy levels into a doublet at ±Em- The suppres- 
sion of local Andreev reflection happens when the width 
Tm of the levels in the doublet (which is finite because of 
leakage into the leads) and the excitation energy E are 
both <C Em- (The relative magnitude of Tm and E does 
not matter.) 

Our theoretical analysis is particularly timely in view 
of recent advances in the experimental realization of 
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topological insulators in two-dimensional (2D) HgTe 
quantum wells [H[ and 3D BiSb crystals [l9[ . Topologi- 
cal insulators are characterized by an inverted band gap, 
which produces metallic states at the interface with vac- 
uum or any material with a normal (noninverted) band 
gap [2o| | . The metallic states are 2D surface states if the 
insulator is 3D, while if the insulator is 2D the metallic 
states are ID edge states. 

These recent experiments [H, HH used nonsupercon- 
ducting electrodes. A superconducting proximity effect 
between Nb and BiSb was reported in earlier work [2l| . 
so that we expect a search for the predicted Q Majorana 
bound states to be carried out in the near future. Antic- 
ipating these developments, we will identify observable 
consequences of the suppression of local Andreev reflec- 
tion, by calculating the shot noise in a 2D topological 
insulator with a superconducting electrode (Fig. [TJ lower 
panel). A similar calculation can be done for the 3D 
case, and indeed our conclusions arc quite general — as 
we will now demonstrate by showing that the Majorana 
Hamiltonian fTJ) directly implies the suppression of local 
Andreev reflection. 

For this purpose write the unitary scattering matrix 
S(E) in a model-independent form, 

S(E) = 1 + 2mW i (H M -E- i-kWW^)- 1 W, (2) 

with W the matrix that describes the coupling of the 
scatterer (Hamiltonian Hm) to the leads. In our case, 
we have 



W ~ 1 w 2 w* 2 ' Hm 





-iE M 



iE M 




(3) 



The expression for Hm is Eq. |T]) in the basis {<I>i, $2} of 
the two Majorana bound states, while W is the coupling 
matrix in the basis {$ e ,i, ®e,2, ®h,i, ^11,2} of propagating 
electron and hole modes in leads 1 and 2. We have as- 
sumed that lead 1 is coupled only to bound state 1 and 
lead 2 only to bound state 2, and we have also assumed 
that the energy dependence of the coupling amplitudes 
Wi can be neglected. (In the exact calculation given later 
on for a specific model, neither assumption will be made.) 
Without loss of generality we can choose the w^s to be 
purely real numbers by adjusting the phases of the basis 
states in the leads. 

Substitution of Eq. §5§ into Eq. @ gives the electron 
and hole blocks of the scattering matrix, 



S ■ 



„hh 



1 + A A 
A 1 + A 



(4) 



which turn out to depend on a single 2x2 matrix A with 
elements 



(The width Tm introduced earlier equals Ti + T2-) Uni- 
tarity of S is guaranteed by the identity 



A + A ] + 2AA^ = 0. 



(7) 



In the limit of low excitation energies and weak cou- 
pling to the leads, this simplifies to 
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for E,Ti < E M . 



(8) 



The scattering matrix s he = A that describes Andreev 
reflection of an electron into a hole has therefore only 
off-diagonal elements in this limit, so only crossed An- 
dreev reflection remains. More specifically, an electron 
incident in lead 1 is transferred to the other lead 2 ei- 
ther as an electron or as a hole, with equal probabilities 
p = Y1Y2/ E M , The probability for local Andreev reflec- 
tion is smaller than the probability p for crossed Andreev 
reflection by a factor (T 1 /T 2 )(E 2 / E 2 M + T\/E 2 M ) < 1. 

Because the probabilities to transfer to the other lead 
as an electron or as a hole are the same, crossed Andreev 
reflection cannot be detected in the time averaged cur- 
rent Ii in lead i, but requires measurement of the current 
fluctuations SIi(t) — liit) — Ii- We consider the case 
that both leads are biased equally at voltage V, while 
the superconductor is grounded. At low temperatures 
T <C eV/ks the current fluctuations are dominated by 
shot noise. In the regime p <C 1 of interest, this noise con- 
sists of independent current pulses with Poisson statistics 
[22| . The Fano factor (ratio of noise power and mean cur- 
rent) measures the charge transferred in a current pulse. 

The total (zero frequency) noise power P = Py, 
with 



dt8Ii(0)8Ij(t), 



(9) 



has Fano factor F — P/el (with I = Y^i^i) ec l ua l to 2 
rather than equal to 1 because the superconductor can 
only absorb electrons in pairs (23|. As we will now show, 
the suppression of local Andreev reflection by the pair of 
Majorana bound states produces a characteristic signa- 
ture in the individual noise correlators Pij . 

The general expressions for J,; and Pij in terms of the 
scattering matrix elements are [24J 



eV 



hh\ 



Ii = T dEll-Kg+n 
h Jo 

2 f eV 

Pij = ^J q dEVij{E), 



(10) 
(11) 



with the definitions 



. _ x (%Cx{E + iT 2 ) -Ems/T^X , . 

\ EmVT^ iT^E + iT,)) ■ W 

We have abbreviated 

Z = E M -(E + iT 1 )(E + iT 2 ), T l = 2iTwl (6) 



(12) 

n x i V(E)=J2^(E)[s v j e k( E )}^ x,V£{e,h}. (13) 
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Substitution of the special form ([4]) of S for the pair of 
Majorana bound states, results in 

Ii = ^r dE(AA%, (14) 
2e 2 r eV 

Pij = eliSij + — J dE [ \Aij + {AA%\ 2 



-\{Arf) 



ij i ] ) 



(15) 



where we have used the identity ([7]). 

We now take the low energy and weak coupling limit, 
where A becomes the off-diagonal matrix ©. Then we 
obtain the remarkably simple result 



el 

Pi, = eh = eh = — , for eV, r 4 < P M . 



(16) 



The total noise power P = V . . P i7 - = 2el has Fano fac- 
tor two, as it should be for transfer of Cooper pairs into 
a superconductor [23[, but the noise power of the sepa- 
rate leads has unit Fano factor: Fi = Pu/eh = 1. Be- 
cause local Andreev reflection is suppressed, the current 
pulses in a single lead transfer charge e rather than 2e 
into the superconductor. The positive cross-correlation 
of the current pulses in the two leads ensures that the to- 
tal transferred charge is 2e. This "splitting" of a Cooper 
pair is a highly characteristic signature of a Majorana 
qubit, reminiscent of the h/e (instead of h/2e) flux peri- 
odicity of the Josephson effect p], Hf| Hsj ■ 

Notice that for any stochastic process the cross- 
correlator is bounded by the auto-correlator, 



\Pj3\ < HPu+Pz 



(17) 



The positive cross-correlation (| 16[) is therefore maximally 
large. This is a special property of the low energy, weak 
coupling limit. There is no inconsistency with the con- 
clusion of Bolech and Demler [llj , that the currents into 
two Majorana bound states fluctuate independently, be- 
cause that conclusion applies to the regime eV ^> Em- 
The duration h/eV of the current pulses is then shorter 
than the time Ti/Em needed to transfer charge between 
the bound states, so no cross-correlations can develop. In 
this high-voltage regime the two Majorana bound states 
behave as independent Andreev resonances, for which the 
noise correlators are known 12711 . 



Pa = eh, /' 
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0, for eV^>E M ,Ti 



(18) 



While the Fano factors of the individual leads Fi = 1 
remain the same, the total noise power P = £^ . P^ = el 
has Fano factor F = 1 rather than F = 2 when the cross- 
correlator Pi 2 vanishes in the high-voltage regime. 

As a specific model that can be solved exactly and 
is experimentally relevant, we consider a 2D topolog- 
ical insulator contacted at the edge by one supercon- 
ducting electrode in between a pair of magnets (Fig. 
[U bottom panel). As discoverd by Fu and Kane Q, 
a Majorana bound state appears at the intersection of 




FIG. 2: Data points: Auto-correlator Pn (circles) and cross- 
correlator Pi 2 (diamonds) of the current fluctuations for the 
model Hamiltonian (|19p . The parameters chosen are Ef = 
0, <f> = 0, mo/Ao = 1, l = 2.3^0, h = h = 3£o. The 
correlators are normalized by eh, to demonstrate the low- 
and high- voltage limits (|16[1 and (|18|) . The dashed and solid 
curves result from the model-independent scattering matrix 
@, with the parameters given by Eqs. PD) and (pH) . The 
dotted curve is the corresponding result for the total noise 
power P = ■ Pij i normalized by e I — e Ii. 



the magnet-superconductor interface with the edge of 
the insulator. The four-component wave function ^> — 
(^ef; ^ej.! ^hfi ®hl) of the edge state satisfies 0] 



rn ■ ct 



- vpa z 
A* 



E F A 
rn a - vpa z 



E, 



* = EV. 



. (19) 

Here p = —ihd/dx is the momentum operator, Ep the 
Fermi energy, v the Fermi velocity, A the supercon- 
ducting pair potential, m the magnetization vector, and 
c = (fx; Gy, &z) the vector of Pauli matrices (acting in 
the space of right and left movers |, J.). 

We set A(x) = everywhere except A = Ao for 
< x < lo- We also set m(x) = everywhere ex- 
cept m — (mo, 0,0) for —l\ < x < and m = 
(mo cos 4>, mo sinci, 0) for Iq < x < Iq + h- We assume 
that | mo | > \Ep\i so that the Fermi level lies in a gap in 
the magnets as well as in the superconductor. The de- 
cay length in the superconductor is the coherence length 
£o = fr^/Ao, while the decay length in the magnets is 
given by A = hv(ml — Ep)^ 1 / 2 . For Ao < £o the only 
bound state at the magnet-superconductor interface is 
the zero-energy Majorana state. 

We have calculated the scattering states for this model 
by matching the 'J's at the opposite sides of the four in- 
terfaces x — —Ii,0,Iq,Iq + h- The resulting scattering 
matrix is then substituted in the general expressions ()10l 
I13[) to obtain the zero-temperature, zero-frequency noise 
correlators as a function of the applied voltage V. Repre- 
sentative results are shown in Fig. [2] (data points) . At low 
voltages we confirm the unit Fano factor and maximal 
cross-correlation of Eq. (Tj"rJ)) , obtained from the model- 
independent scattering matrix ^ . Also the crossover to 



4 



the conventional high- voltage regime (fT8|) of independent 
resonances is clearly visible. 

For a quantitative comparison of the two calculations 
we need the splitting and broadening of the Majorana 
bound states in the tunneling regime 1%, 1% ^> Ao, Iq 3> £o- 
We find 



E 



M 



e -i /io cos 



Epl 

hv 



arctanf 



Ep\ , 
hv ' 



_ e -2h/\ 



r, = e 



(1 - E 2 F /ml) 



2hv 
Co + A 



2hv 

(20) 
(21) 



Notice that the level splitting can be controlled by vary- 
ing the angle <p between the magnetizations at the two 
sides of the superconductor [28[ . In Fig. [2] we use these 
parameters to compare the model-independent calcula- 
tion based on the scattering matrix @ (curves) with the 
results from the model Hamiltonian (p^9|) (data points), 
and find excellent agreement. 

The setup sketched in Fig. Q] might be realized in a 
HgTe quantum well (l8j . The relevant parameters for 
this material are as follows. The gap in the bulk insu- 
lator is of the order of 20 meV and the magnetic gap 
can be as large as 3 meV at a magnetic field of 1 T. 
The smallest energy scale is therefore the gap induced 
by the superconductor, estimated [26[ at A = 0.1 meV. 



With hv = 0.36 meV • /mi this gives a superconducting 
coherence length of £o = 3.6 /jm, comparable to the mag- 
netic penetration length Ao at a field of 0.03 T. For the 
calculation in Fig. [5] we took £o = Ao and then took 
the length Iq of the superconducting contact equal to 
2.3 £o — 8 /mi, and the lengths h,h of the magnets 
both equal to 3£o — lljum. The level splitting is then 
E M = 0.1 A = 10/teV = 100 mK. At a temperature 
of the order of 10 mK we would then have a sufficiently 
broad range of voltages where ksT < eV < Em. 

In conclusion, we have demonstrated the suppression 
of local Andreev reflection by a pair of Majorana bound 
states at low excitation energies. The remaining crossed 
Andreev reflection amounts to the splitting of a Cooper 
pair over the two spatially separated halves of the Ma- 
jorana qubit. This nonlocal scattering process has a 
characteristic signature in the maximal positive cross- 
correlation (P\2 = Pn = P22) of the current fluctua- 
tions. The splitting of a Cooper pair by the Majorana 
qubit produces a pair of excitations in the two leads that 
are maximally entangled in the momentum (rather than 
the spin) degree of freedom, and might be used as "flying 
qubits" in quantum information processing. 
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by the Dutch Science Foundation NWO/FOM. 
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